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We study the phenomenon of decay of a supercurrent in a superconducting thin 
film due to the spontaneous homogeneous nucleation of quantized vortex-antivortex 
pairs in the absence of an external magnetic field and in the presence of a periodic 
pinning potential. We describe the vortex nucleation by means of a Schwinger-type 
path-integral calculation including the effects of quantum dissipation. The over- 
damped case is treated exactly, whilst the case in which a periodic array of pinning 
centres is present is dealt with by means of a renormalization group (RG) approach 
in frequency space. This yields an approximate but very appealing analytic result 
for the dependence of the vortex nucleation rate on the externally-driven supercur- 
rent. In this formula the rate dependence displays oscillations which are connected 
to the pinning periodicity and correspond to the vortex nucleation length probing 
over the pinning sites as the current density is varied. Our RG treatment describes 
the vortex nucleation process in both the confined and the mobile phases of the 
dissipation-driven localization transition characterising the motion of dissipating 
quantum particles in a periodic potential. 



1 Introduction 

There has been a strong revival of interest, recently, in the physics of magnetic 
vortices in type II and high-temperature superconductors □. Most research 
efforts have been devoted to phenomena relating to the nature of the mixed 
phase of a superconductor in some externally applied magnetic field and su- 
percurrent. Issues connected with the pinning of the flux lines by defects have 
been widly studied. 

We □, as well as Ao and ThoulessB and Stephen □, have addressed the 
problem of the quantum dynamics of vortices in the absence of an external field 
but in the presence of an externally driven supercurrent, quantum dissipation 
and pinning. This leads to the decay of a supercurrent, or a residual zero- 
temperature resistance in the superconductor. Whilst most of the dissipation 
seems to be ascribed to vortices tunneling in the sample from the edge, an 
interesting novel possibility also explored by us in depth is that of a residual 
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resistance arising from spontaneous vortex-antivortex pair creation in the bulk 
of a thin film. This is the mesoscopic counterpart of electron-positron pair 
production of two-dimensional (2D) quantum electrodynamics (QED) in the 
presence of static e.m. fields, which in a superconductor arise from the static 
and velocity-dependent components of the Magnus force acting on the vortices. 
Exploiting this analogy with QED, a powerful "relativistic" quantum field 
theory approach has been set up to study vortex nucleation in the 2D geometry 
in the presence of quantum dissipation and of pinning potentials. The central 
result is that the nucleation rate T has a strong exponential dependence on the 
number current density J, given by 

T^ 2 Veff J' 1 eM-VefA/^J 2 } (1) 

Here r] e f / is an effective viscosity coefficient as renormalised by the magnetic- 
like part of the Magnus force, and £qr is the rest- or nucleation-energy of 
a single vortex as renormalized by screened Coulomb interactions and (fake) 
Landau-level corrections. This exponential dependence would make the vortex 
nucleation (folded, e.g., into the sample's resistance) observable in a rather 
narrow range of J-values. Thus, normally the superconductor is essentially 
resistance-free. However, the high values of J that can be reached in the high- 
T c materials make the possibility of observing pair creation in static fields 
within reach for thin films. One particular feature that would uniquely relate 
the residual resistance to the phenomenon of spontaneous vortex-pair creation 
is the presence of oscillations in the J-dependence of T( J) in case a periodic 
pinning potential is artificially created in the film. These oscillations are in 
fact strictly connected to the pinning-lattice spacing d — 2i:/k of the periodic 
potential (we assume a square lattice), e.g. 

2 

U(q(t)) =U J2^- cos (kq a (t))} (2) 

o=l 

acting on the nucleating vortex-pairs described by a coordinate q. 

The problem of quantum dissipation for a particle moving in a periodic 
potential has some interesting features in its own right Wail. It is charac- 
terised by a localization phase transition driven by dissipation; accordingly, 
two phases can occur depending on whether the dissipation coefficient El r] 
is greater (confined phase) or smaller (mobile phase) than a critical value 
Vc = k 2 /2n = 2ir/d 2 . This localization transition is described by a Kosterlitz- 
type renormalization group (RG) approach, yet with some important differ- 
ences that will be recalled below. We have implemented the RG approach for 
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the evaluation of the dependence of the spontaneous nucleation rate of vortex- 
antivortex pairs on the external parameters for our own quantum dynami- 
cal system. A remnant of the dissipation-driven phase transition is observed 
and the pair production rate T can be derived in both phases by means of 
a frequency-space RG procedure leading to observable current-oscillations if 
t] > Vc- 



2 RG approach to dissipative localization transition 

First, we briefly recall the RG description of the localization transition driven 
by quantum dissipation 0. The effective action for a particle diffusing in a 
periodic potential and subject to quantum dissipation of the Caldeira-Leggett 
typeta is, in Fourier frequency space: 



S= £(q) =T^{-mujl + -?]\uj n \}q a (uj n )q a (-uj n ) + dtU(q) 
Jo „ 1 1 Jo 



(3) 



where m is the mass of the quantum particle and r\ the phenomenological fric- 
tion coefficient. In the low-frequency limit the effects of inertia can be neglected 
and the problem would acquire the same phenomenology as for the sine-Gordon 
model (in (0+l)-dimensions), except for the peculiar aj„-dependence of the 
propagator reflecting the broken time-reversal symmetry of quantum dissipa- 
tion. When the RG procedure is applied to Eq. (||) a renormalization of the 
potential amplitude Uq occurs, but not of the friction coefficient 7/ since only 
local operators in the time variable can be generated within a RG transforma- 
tion. In terms of the dimensionless parameters (f2 is a large frequency-cutoff) 
U = Uq/VL and a = 2irri/k' 2 , the RG recursion relations read 

with e~ e the frequency-scale renormalization parameter. These have the simple 
solution 

U(£) = U(0)e {1 ~^ Me , a{£) = a(0) (5) 

displaying the localization transition for 77 = r\ c = k 2 /2n = 2ir/d 2 . The poten- 
tial's amplitude vanishes under a RG change of time scale for r\ < r) c , but for 
77 > it tends to diverge and the RG procedure must be interrupted. Unlike 
in the Kosterlitz RG scheme, this cannot be done unequivocally in the present 
situation, for there is no true characteristic correlation time or frequency owing 
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to the fact that one never moves away from the neighbourhood of the critical 
point r) c . An alternative strategy for the confined phase is to resort to a vari- 
ational treatment 13, which dynamically generates a correlation time. In this 
procedure the action of Eq. (||) is replaced by a trial Gaussian form (neglecting 
inertia) 



where M 2 is determined by minimising self-consistently the free energy F tr + 
(S — S tr )tr- This leads to the equation 

having a solution M 2 ^0 only in the confined phase (77 > rj c ), since introducing 
the cutoff f2 in the (continuous) sum over frequency modes w„ = 27rn/r, the 
equation for M 2 leads to (for M 2 — » 0) 

M 2 = ^^|J =r,Q» (8) 

This spontaneously generated "mass" interrupts the divergent renormalization 
of the periodic potential amplitude Uq, which in the RG limit £—>oo tends to 

fY""' „ /V+i/»* w " 



Here, we have put fl = 2irn* /t = n*U\ and /1 = M 2 /Clrj. 



3 RG treatment of vortex-antivortex pair-creation in the presence 
of a periodic pinning potential 

We begin by recalling the need for a relativistic description of the process. This 
leadscl to a Schwinger-type formula for the decay of the "vacuum" , represented 
by a thin superconducting film in which static e.m.-like fields arise when a 
supercurrent is switched on. The quantum fluctuations of these fields are 
vortex-antivortex pairs, nucleating at a rate given by 

T 1 poo j p pt 

T 2 = L^ Im I T e ~ £ °J p 9(*) ex P{- J dtc ^ ( 10 ) 
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where L 2 T is the space-time volume of the sample and the vortex-nucleation 
energy (suitably renormalised by vortex-screening effects). Also 



AA^ - \nnFiJ.vqv + V(q) 




(11) 



is the Euclidean single-particle relativistic Lagrangian, incorporatingJhe pin- 
ning potential V(q) — 2£ U(q) and the Caldeira-Leggett mechanismEl. In the 
absence of the pinning potential, the relativistic action is quadratic and the 



path integral in Eq. (10) can be evaluated exactly. The leading term in the 
expression for T follows from the lowest pole in the r-integral and this can 
be obtained exactly in the (non-relativistic) limit in which mi = = ^— > 0, 
with - — £ /m^oo playing the role of the square of the speed of light. The 
result! is Eq. (0). 

We now come to the evaluation of T in the presence of the periodic po- 
tential, which calls for the RG approach of Section 2. Integrating out the 
Euclidean "time" -like component (73 (t), we reach a formulation in which the 
electric-like and the magnetic-like Magnus field components are disentangled. 
In terms of Fourier components, dropping the magnetic-like part and for 7— >0: 



f dt£ E {q) =rV{^K| -E 2 8 al }q a {u n )q a (-Lu n )+ f dtV(q) (12) 

with E = 2tlJ the electric-like field due to the supercurrent donsity J. Wc 
have shown □ that the only role of the magnetic-like field is to renormalizc 
the nucleation energy and the friction coefficient, hence our problem amounts 
to an effective one-dimensional system in the presence of E and dissipation. 
The evaluation of the Feynman Path Integral (FPI) proceeds by means of 
integrating out the zero- mode, q~o, as well as the high-frequency modes q n with 
n > 1, since again the leading term for T in Eq. ( |l0|) comes from the divergence 
of the FPI associated with the lowest mode coupling to E. The effect of q n 
with n > 1 is taken into account through the frequency-shell RG method of 
Section 2, leading to a renormalization of the amplitude Vo = 2£qUq of the 
(relativistic) pinning potential. The renormalization has to be carried out 
from the outer shell of radius Q to lu± = 2tt/t. In the mobile phase (r/ < rj c ) 
this implies e l — £It/2tt = n* in Eq. (||), with (from the leading pole of the 
FPI) r = irrj/E 2 (— >oo for relatively weak currents). In the more interesting 
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confined phase (77 > ?y c ) we must integrate out the massive n > 1 modes 
with a Lagrangian £(q n ) = t [-^r]\uj n \ + \M 2 — E 2 ) g n g*. This leads to an 
additional, entropy-like renormalization of the activation energy £oi?., beside 
the renormalization of Vo- We are therefore left with the integration over the 
modes q and q~\ , with a renormalised potential 




dtV R (q ,qi(t)) = V t - Vqr / dtcos(k(q + qi (t))) 
Jo 
^V t - V OR TJo{2k\ qi \) cos(fego) (13) 

Here, Jq is the Bessel function and the renormalised amplitude Vqr is 



In Eq. (|l^) the phase of the <?i mode has been integrated out, allowing us to 
integrate out the go-mode exactly; this leads to the expression 



p />oo rOO 

^2= Im J dTM(T)e~ {£ ^ +Vo)T J dW^e-^-z 2 ^ 2 I Q (V 0R r J (2k\qi\)) 

(15) 

where Iq is the modified Bessel function. It is clear that the singularity from 
the gi-integral occurs at r = irrj/E 2 ; evaluating the normalization factor A/"(r), 
we finally arrive at 



r = r Q K(j) (16) 

^)=e(l + M)(l + ^)/o(gfjo(2^) 

where T is given by Eq. ([l]), there is a further renormalization £qr^£qr + 
\q and we have set E = 2irJ. is a nucleation length, which is in first 
approximation given by 



l N- 



4tt 2 J 2 4tt 2 J 2 



TV J 



(17) 



and corresponds physically to the distance a vortex and antivortex pair must 
travel to acquire the nucleation energy £qr. The presence of the Jq(2Mn) 
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argument in the correction factor K{ J) due to the pinning lattice thus gives rise 
to oscillations in T(J) (hence in the sample's resistance) through the parameter 
2fc£jv = <±kIn jd. Vortex nucleation is therefore sensitive to the corrugation of 
the pinning substrate. However, these oscillations should be observable only 
in the confined phase, r\ > r) c , where interrupted-renormalization prevents the 
prefactor in front of the Jo{x) oscillating function from becoming too small for 
relatively small current densities. 
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